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Abstract 

The generalization of the Glauber formula for cumulative production events is derived. On 
its basis the multiplicity distribution in such events is related to the one in the minimum 
bias events. As compared to the rare events of type C, the formula involves a shift in the 
arguments determined by the multiplicity from a collision with a cluster of several nucleons. 
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1 Introduction. 

In [1,2] a simple formula was derived which relates the multiplicity distribution Pc{n) in a 
rare event of the so-called type C to the minimum bias (MB) multiplicity distribution P{n): 

Remarkably the right-hand side of (1) does not depend on a particular nature of the rare 
process. So (1) gives a universal probability for rare events, independent of the nature of the 
trigger process. This universal formula was tested in several different rare processes for hh 
and hA collisions like W and Drell-Yan dilepton pair production. 

In this paper we would like to discuss the possibility of applying (1) to cumulative events in 
hA collisions. In deriving (1) it was assumed that the geometric properties of the hA collision 
in the rare event were the same as for MB events. Prom this point of view production 
of cumulative particles, that is, in the kinematical region outside the one allowed for the 
scattering on a nuclcon at rest, is different, since it involves collisions with two or more 
nucleons at once. Thus it is not evident that (1) is applicable to cumulative events, although 
they are certainly rare, experimentally. 

Our results show that under certain reasonable approximations one obtains an expression 
for the multiplicity distribution in a cumulative event quite similar to (1), except for a certain 
shift in the arguments. The shift is due to a possible change of the mutiplicity in a collision 
with a cluster of several nucleons. In the colour string approach, clustering of nucleons leads 
to fusion of strings and their percolation, which tends to diminish the nuiltplicity [3,4] The 
mentioned shift in the arguments can be observed experimentally and thus can serve as a 
signature of string fusion, once the multiplicity distributions in cumulative and MB events 
are compared. 



2 The Gribov-Glauber formula for cumulative events 

To study the multiplicity distribution in a cumulative event, one needs a suitable general- 
ization of the standard Glauber formula for hA collisions. The derivation of the standard 
Glauber hA amplitude for high energies made by V.Gribov [5], expresses the hA amplitude 
via hN amplitudes with the nuclcon at rest. The corresponding kinematical region then ev- 
idently coincides with the one for the scattering on a nucleon at rest. Cumulative events, 
in contrast, involve collisions with several nucleons at once. Such multinucleon collisions 
were considered earlier in relation to the loop contribution in nuclear-nuclear collisions [6] or 
string fusion [7]. It is remarkable, however, that one can naturally introduce them strictly 
following the Gribov approach, provided one takes into account a more general structure for 
the high-energy part of the hA amplitude. Here we only present the main steps and the final 
formula. One can find some details in the Appendix. 

The basic novelty in deriving the Glauber amplitude for cumulative events is taking into 
account that the high energy part H of the hA scattering amplitude (see Figure) , apart from 
the usual terms which depend on all momenta transferred to the nucleus qi separately, may 
include other terms which depend on their partial sums, say, qi + q2, ?3 + ?4 + Qb etc. As we 
shall see, such terms describe simultaneous interaction of the projectile with clusters of two, 
three etc. nucleons. Correspondingly, we shall call these sums of the transferred momenta 
clusters for brevity. In the general case, clustering implies a subdivision of all transferred 
momenta into various groups (clusters) of one, two or more momenta. The high- 

energy part can be split into terms which depend only on the total momentum of each cluster. 
Some of the clusters may contain equal number of nucleons, say, q\ + q2 and q^ + g4 with 
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two nucleons each. The order of such clusters has no influence on the resulting contribution 
to the amplitude. For this reason it is convenient to characterize the clustering structure 
by occupation numbers v^i k = 1, 2, .... telling the number of clusters containing k nucleons 
(that is, k transferred momenta qj). Evidently for the n-fold interaction with the nucleons 
one should have 

^kuk = n (2) 

k 

and, of course, all with k > n are zero. The standard situation when the high-energy part 
depends on all qj separately corresponds to z^i = n and ffc = for A: > 1. 

Taking the possibility of clustering into account, we present the high- energy part as a 
sum of contribution from all possible clusters: 

H{q^,...qn)=XlY.Su,H,,{qi,...qn) (3) 

k I'k 

Here Hnj^i^qi, ■■■qn) = H^j^{ki, K2, ■■■■Ku) is supposed to depend only on the total cluster mo- 
menta Kj, j = 1, 2, ....1/ where 

= ^^k 
k 

is the total number of clusters. We also assume that the initial momenta qj arc distributed 
among the clusters in some fixed way, say, in the order of the growing cluster size. Various 
redestributions are then taken into account by the evident symmetry factor 

Eqs. (3) and (4) present the essence of our generalization of the Gribov derivation. Then, 
following the standard procedure of integrating over the longitudinal momentum transfers by 
means of closing the contour around the right-hand singularities of the high-energy part and 
retaining the lowest mass singularities, one readily finds a contribution to the hA amplitude 
from a given clustering vj^ 

iA., = S,, f d\Jl[iakn{r)Y^ (5) 

k 

Here are the scattering amplitudes for a collision of the projectile with a cluster of k 
nucleons at rest. Their normalization is chosen to include an extra factor (2s)~^(2m)^~''^ 
as compared to the standard relativistic one, where m is the nucleon mass and s the cm 
energy squared. With this normalization for fc = 1 (and for the total hA amplitude) twice 
the imaginary part directly gives the total cross-section. The generalized profile functions 
are defined as 

Tfc(r) = / dzp''{r,z)dz (6) 



where p(r, z) is the 3-dimensional nuclear density. As expected, they give the probability to 
find k nucleons at the same point in the transverse plane. 

To find the total hA amplitude we have to sum (5) over all cluster structures for fixed n 
and then sum over all n with the standard factor C^. This reduces to summing over all 
with the only restriction 

Y,kyk<A (7) 

k 

In this way we obtain the hA amplitude in the form 

A = t d'^rA{r) (8) 
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where the hA amplitude for fixed impact parameter is given by 



and the summations are restricted by condition (7). 

The standard Glauber formula is obtained if all = for A; > 1. In the general case we 
cannot do the summations explicitly for arbitrary A. However they are easily done for a very 
heavy nucleus with A» 1. In this limit we can drop condition (7) and approximately take 

~ TT A'"''' 

iA-j:kk^ky.-\^ 

Then the summation over each leads to an exponential factor and we obtain 

iA{r) = exp ^ i^k^-^^ - 1 (10) 

The subtracted unity corresponds to a term with all z^^ = and thus without any interaction. 
Evidently, neglecting all amplitudes with k > 1 one obtains the standard Glauber expression 
for a heavy nucleus. 

Note that ~ ^^/^-^ and so A''Tk ~ A^/^. As a result all terms in the exponent in 
(10) have the same order A^^^ independent of the number of nucleons in the cluster k. This 
means that the relative weight of the multinucleon interactions is essentially independent of 
A. Their smallness is determined not by the large nuclear volume as a whole but by the large 
internucleon distance Rq in the nucleus as compared to the strong interaction radius Rg, 
which determines the mag nitude of a^. Indeed A'^Tk ~ A^s^i-Sfc^ relative weight 

of the interaction with k nucleons at once is of the order of the dimensionless parameter 
akRl~^^ ~ {Rs/Rq)^'^~^- Taking Rs ~ l/rup ~ 0.27/ one finds that raising the number of 
nucleons in the cluster by unity leads to a damping factor of the order (1/4)^ = 1/64 

To calculate the contribution of multinucleon interactions one has to know the multinu- 
cleon amplitudes a^. They only weakly depend on the energy ,s. Although being on-mass-shell 
quantities, they do not seem to be directly measurable. An obvious way is to determine them 
from simpler processes, like the scattering on the deuteron, triton etc, and then use the found 
values for studying the cumulative kinematical region in interactions with other nuclei. 

3 The AGK rules, cross-sections and multiplicities 

To simplify our notations we define tjt(r) = Tk{r)/k\ Then at fixed r 

•^«^S(i^&OT^nKMoi- (11) 



The twice imaginary part of this amplitude gives the total cross-section on the nucleus at 
fixed r 

a^ir) = -iA - {lAT = Y: ^ f , ^Urr[-{iakr ' (12) 

To obtain the AGK rules we have to split this into terms corresponding to a given number 
of inelastic collision with the nuclear clusters. We use a trivial identity 

iak + (iak)* + 2Imafc = 
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It demonstrates that for a given discontinuity of the diagram (cutting of the diagram) each 
interaction may be either cut or taken uncut on both sides of the cut. We present in (12) 

- {iakT'' - {iak^ = (iak + (ia^)* + 2Imafc)'^'= - {iakY" - {iaky"" (13) 

Terms proportional to (2Ima)fc)^'= on the right-hand side correspond to ^^-fold inelastic col- 
lision with clusters k. They are 

C^^H2Imafe)«H^afc + {iakTr'^' = Clla"^' {-1)"'-^" 

where we denoted 

(jjfc = 2Im Ofe 

Terms without 2Im represent the diffractive part of the cross-section 

Summing over all sets of clusters we get the cross-section at a given r corresponding to 
^jt inelastic interactions with cluster k in the form 

^ ? (J^E;OT^ncft(.AM)"'(-i)-^> (14) 

where of course z^fc > ^fc. Changing the summation variable we rewrite this as 

(15) 

This formula expresses the AGK rules in the general case. It can easily be simplified for a 
heavy nucleus when in the same way as in deriving (10) we easily obtain 



a 



A 



(^) = il Tl exp -2^afc^tfc(r-) (16) 

k \ k J 



in an obvious generalization of the well-known formula for collisions with separate nucleons 
inside the nucleus. 

The multiplicity observed in an inelastic collision with a cluster of k nucleons is not 
known apriori. Its magnitude depends on the dynamics of the interaction with the cluster. 
In a model in which this interaction is realized via creation and decay of colour strings, the 
multiplcity fik depends on the interaction between the strings. With no interaction, the 
nucleons, although located at the same longitudinal point, produce strings independently so 
that = kfii. In the opposite case all strings fuse into a pair and //^ = /xi. 

To find the multiplicity in the event with inelastic interactions with cluster k one has 
first to find the corresponding inclusive cross-section. This corresponds to substituting in 
(15) one of the cross-sections by /UfcCTfc and then summing over all such substitutions. The 
multiplicity is obtained by dividing the result by the cross-section (15) itself. This gives a 
simple formula 

li^^"'^ = Y.^k^ik (17) 

k 
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4 Rare events. 

As discussed, the probability for the nucleons to form a cluster is rather small. So in the 
cumulative event, in which at least one cluster with a given fc > 1 interacts with the projectile, 
in the first approximation one can neglect all configurations with two or more clusters with 
A; > 1. We shall assume that the nucleus has a sharp edge and a constant density inside, so 
that 

p{r,z) = {l/VA)e{R\-r^ -z^) (18) 
where Ra = A^^^Rq is the nuclear radius and Va = AVq is the nuclear volume. Then we find 

with T(r) = Ti(r) the standard nuclar profile function. 

With these simplifications and neglecting the contribution of clusters in the sum over Uk 
in (15) we find the cross-section for a single inelastic interaction with a cluster k and m — 1 
interactions with separate nucleons at a given impact parameter r as 

(7$"''=^ = 71 r^^^^TTT —^^{aT{r)r{l - aT(r))^-'"-'=+i (20) 

^ {A - m - k + l)\{m - \)\ k\vX~^ 

For a heavy nucleus this transforms into 

(m,fe) _ 1 CTfc 



{m-l)\k\V^ 



l^{aAT{r)re-^^^^''^ (21) 



Let us compare this cross-section with the one for m inelastic interactions with isolated 
nucleons (without clustering) 

"S"^ = (A-1l)!m! ("^("^^"^^^ - -^(^))^"'" (22) 
which for a heavy nucleus goes into 

aS") = _l(^^r(r))'"e-^^^('^) (23) 

The comparison of (21) and (23) gives a relation 

a^^'"^ = a,ma^ (24) 

where we have defined 

Eq.(24) means that the probability Pkim) to find apart from a single interacting cluster k 

also m — 1 interacting separate nucleons, m = 1, 2, is related to the probability Pirn) to 
find m interacting separate nucleons and no clusters at all according to 

Pk{m) = CkakmP{m) (26) 

Since the total probability to have any number of separate interactions with the nucleons 
apart from an interaction with the cluster is unity, the sum of (26) over all m should give 
unity, which determines the constant c^. We then obtain 

PJrn) - ""'^(""'^ - "^-^(""^ (27) 



7 



This relation has the same form (1) as obtained in [1,2] for rare events of the so-called C-type. 
It is remarkable that the relation is independent of the cluster characteristics {k and ak), the 
only requirement being that the interaction with it should be rare. 

However the resulting observable consequences will be different. The point is that passing 
to observable multiplicities, we have to introduce the multiplicity for the interacting cluster 
Then the total multiplicity M, observed in the event in which, apart from the cluster, 
m — 1 isolated nucleons interact will be 



M = + {m — l)n = m/j, + — ^ = m/j, + (28) 

where /i = ni is the multiplicity for the interaction with a isolated nucleon. So, in terms of 
the multiplicities, the relation (27) implies 

, , (M- Afc)P(M- Afe) 

This is our final result. It shows that, applied to rare cluster interactions, the approach 
of [1,2] leads to a very similar relation, in which however the argument is shifted on the 
right-hand side, the shift depending on the multiplicity for the given cluster. As mentioned, 
the latter depends on the assumed dynamics of colour strings in the cluster. Without any 
interaction between them — ^ = [k — 1)^ and the shift is maximal. In the opposite case 
when the cluster behaves in the same manner as a single nucleon 11^ = 11 and there is no shift 
at all. Thus experimental studies of the relation (29) can shed light on the strength of the 
interaction between colour strings and the probability of their fusion and percolation. 

In terms of averages for a cumulative event involving a cluster with k nucleons and () 
for MB events, one gets from (29) for the average multiplicity 

{M)k = J^ + ^k (30) 

and for its dispersion squared 

Dk -{M)k- (M), - - (31) 

One observes that the dispersion of the multiplicity does not depend on the the cluster 
properties and is the same as for rare events of the type C found in [1,2]. As discussed 
there it is smaller than for MB events, so that production of cumulative particles triggers 
sharpening of the multiplicity distribution. On the other hand, the average multiplicity in the 
cumulative event generally results still greater than for rare events of the type C (given by the 
first term in (30)) due to the extra term A/.. Only in the limitiung case when A^ = 0, which 
in the string language corresponds to a very strong fusion, the two multiplicities coincide, 
remaining greater than for MB events. 



5 Conclusions 

The standard Glauber- Gribov derivation of the hA-amplitudc in terms of hN amplitudes 
has been naturally extended to include cumulative effects. In this way the hA amplitude is 
obtained in terms of the amplitudes for the scattering of the projectile hadron off clusters of 
k nucleons, k = 1, 2, .... The resulting formula has the same structure as the original Glauber 
formula. By means of the AGK cutting rules cross-sections for a given number of inelastic 
collisions with clusters of k nucleons are found. Assuming a collision with a cluster of more 
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than one nucleon to be a rare event and thus neglecting contributions from two and more 
colhsons with such clusters, the assocated multiplicity distribution is expressed in terms of 
the MB multiplicity distribution. The obtained formula has the same universal form as (1), 
except for a shift in the arguments. The shift is determined by the difference between 
the multiplicities in the scattering with a cluster of k nucleons and with a single nucleon. 
Observation of this shift potentially gives a possibility to study the multiplicity coming from a 
cluster of nucleons. Although this does not seem to be easy experimentally, the shift could be 
studied in the forthcoming experiments at RHIC and LHC, using, for instance, the BRAHMS 
detector. 
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7 Appendix. The hN amlitude with multinucleon interactions 

The hA amplitude with several (n) rescatterings can be represented by a generic diagram 
shown in Figure. It divides into a part related to the structure of the nucleus and the high- 
energy part represented by the blob H. Separation of the nuclear part is, in fact, standard 
and follows the original approach of V.Gribov [5]. We briefly repeat it for self-consistency 
and to introduce the necessary notations. The latter are as follows. We denote the nucleus 
momentum as Ap and take the nucleus at rest: p = 0. The nucleon momenta before (after) 
the interaction are denoted ki {k'^ = ki + qi). The spectators correspond to i = n + 1, ...A. 
For them ki = k[ and qi = 0. Also ^(/j = 0. The projectile momentum is denoted I. 
The expression for the amplitude corresponding to the diagram in the Figure is 

^-4= / n^7^m)^(^;) n ^P{k,)im)ink'i)iH{l,kj,k'^) (32) 
j=2 y^^r (27r)^ .^^^^ (27r)^ 

In this expression H is the above mentioned high-energy part and P{k) is the propagator of 
the nucleon with momentum k: 

= m2 - fc2 _ ^0 - a2 + k2 _ 2mko - iO ^^^^ 

where m is the nucleon mass and Aa'^/{2m) is the nucleus binding energy. The vertex T 
describes the transition of the nucleus into A nucleons. We have chosen the momenta of the 
first nucleon ki and k[ as dependent variables, although in future we shall use a different 
choice also. Evidently ki = — J2t=2 ki and similarly k'^ = — J2t=2 K- 

Standardly one starts by the integration over the zero components of the momenta. Since 
the poles coming from the propagators of the first nucleon all lie in the upper half-plane, we 
can integrate over kio or /c'q, i = 2, ...A, just taking the residue at the pole of the corresponding 
propagator P{ki) or P{k[). Each such integration provides a factor 27r/(2m). The two 
propagators of the active nucleon in the initial and final state together with the factors 
iV{ki)iV{k[) combine into a product of two nuclear wave functions 

so that the expression for the amplitude (32) becomes 
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The next step is to pass to the coordinate space to simphfy the dependence on the wave 
functions. We present 

<^(^^) = / n (^^(^^) ^^p(-^ E '^j^^) (35) 

and then integrate over the transverse momenta. The high-energy part H does not depend 
on the transverse momenta of the nucleons in the high-energy hmit. So the integration over 
them is trivial, since all the relevant dependence is concentrated in the exponentials. We 
readily obtain 

jA- f TT TT ^'^^^ 

iH{lz,kjz, Qjz) n exp(-z ^ kzj{zj - z'^)) ]J dz[ exp(i ^ qzj{z'j - z[)) 

A-l 

d^ri ]J d'^rjipin =r2 = r3 + ...rn, rj; Zi)il){ri = r2 = rs -h ...r„, ry, z[) (36) 

j=n+l 

As expected on the physical grounds, the n nucleons which take part in the rescattering have 
to be taken at the same impact parameter. 

Now we have to integrate over the longitudinal momenta. Evidently the integrand docs 
not depend on the kz of the spectators. So these interactions are done trivially and add a 
factor 

A-l 

n 2T,5{zi-z[) 

i=n+\ 

Apart from the factor {2it)^~'^~^, these 5 functions convert the double integration over z 
into a single one for the spectators. Together with the integration over their transverse 
coordinates this turns the product of the wave functions into the nuclear p-matrix for n -|- 1 
nucleons taking part in the interaction. 

The high-energy part H does not depend on kjz but only on qjz- This property follows 
from the relativistic invariance, due to which H depends on products like Ikj and Iqj. With 
small spatial componets of kj, the products of the first tipe all arc approximately equal to Ip 
and only the second are real variables. So one can also trivially integrate over kjz-i j — 1, ...ti^ 
which gives a factor (27r)" and puts zj = Zj, j = 1, ...n We get 



= / n l^^^^^^iH{lz,qjz)Y[dZieyiY>{i'Yj(lzj{zj - zi))d^rp{^ 

(37) 

At this point all nuclear effects are taken into account by the nuclear p-matrix taken at 
interaction points (r, Zj ) . The final task is to integrate over the transferred longitudinal 
momenta qjz- To do this we take into account the generalized structure of the high-energy 
part as a function of the transferred momenta, discussed in Sec. 3 and presented in Eqs (3) 
and (4). 

Let us find the contribution to the hA amplitude with n interactions (7) from a term in 
(3) with given v^. All integrations over qj can be divided into v — 1 integrations over the 
total cluster momenta K2i---Ku (^i = — X]j=2'^j) ^^^^ ^ji without, say, the first q 

from each cluster. The integration over the latters is trivial, since, the high-energy part does 
not depend on them. It will make all longitudinal points Zj equal within each cluster, in 
accordance with the intuitive definition of a cluster or simultaneous interaction with several 
nucleons. 
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To write down the obtained expression in an understandable way, at this point we make 
the usual simplifying assumption about the structure of the nuclear p matrix: that of absence 
of correlations and consequently of factorization 

p{rzi,rZ2, ....rzn\rzi,rz2, ...rzn) = Wp{rzj) (38) 

where p{rz) is the usual nuclear density (normalized to unity). Then the contribution to the 
amplitude from will be given by 

iA = 5.,(2m)-"+i / n -^iH^d^j) n dzi exp(i^ «,(z,- - z,))d'r [] p'^{rzj) (39) 

j=2 i=l j=l 

where kj is the number of nucleons in the cluster at point zj 

This expression can now be treated in the same way as in the original Gribov approach. 
Instead of Ki we introduce u — 1 cumulative longitudinal momentum transfers ti — I X^fci ^j- 
This gives a factor {2iy~'^ . The Feynman integration contour over each of ti can be closed 
around the right-hand side singularities. In the spirit of the Glauber approximation, from 
these singularities we retain one the pole singularity at = corresponding to the single 
nucleon intermediate state. Taking the residue at this point makes the exponential factor 
in (12) equal to unity. The residue of iH^^^ itself is a product of v on-mass-shall forward 
connected amplitudes ia^ for the interaction of the projectile with k nucleons at rest, where 
k is the number of nucleons in a cluster. After that the final integration over the cluster points 
Zj becomes trivial and leads to a standard profile function and its evident generalization. We 
obtain 

iA,, = S,,{4ml)'-''{2mr-^ ( <fr \{[iakn{r)r^ (40) 

k 

where 

Tfc(r) = j dzp''{r,z)dz (41) 

is an obvious generalization of the profile function. Note that Ami = s where s is the standard 
cm. energy squared. It is convenient to pass to reduced amplitudes 

Qfc ^ o /o°\a-i ' dimojk = l-3fc (42) 
2s(2m)'^ 

and similarly for A. After that one obtains Eq. (5). 
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9 Figure caption 

The hA scattering amplitude. Double lines show the nucleus target. 



